Non-Drude universal scaling laws for the optical response of local Fermi liquids 
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We investigate the frequency and temperature dependence of the low-energy electron dynamics in a 
Landau Fermi-liquid with a local self-energy. We show that the frequency and temperature dependence of 
the optical conductivity obeys a universal scaling form, for which an explicit analytical expression is obtained. 
For the optical conductivity and the associated memory function, we obtain a number of surprising features 
which differ qualitatively from the Drude-model and are universal characteristics of a Fermi liquid. Different 
physical regimes of scaling are identified, with marked non-Drude features in the regime where ha> ~ fc B 7\ 
These analytical results for the optical conductivity are compared to numerical calculations for the doped 
Hubbard model within dynamical mean-field theory. For the 'universal' low energy electrodynamics, we 
obtain perfect agreement between numerical calculations and analytical scaling laws. Both results show 
that the optical conductivity displays a non-Drude 'foot', which could be easily mistaken as a signature of 
breakdown of the Fermi liquid, while it actually is a striking signature of its applicability. The aforementioned 
scaling laws provide a quantitative tool for the experimental identification and analysis of the Fermi-liquid 
state using optical spectroscopy, and a powerful method for the identification of alternative states of matter, 
when applicable. 



PACS numbers: 78.20.Bh, 78.47.db, 72.15.Lh 

I. INTRODUCTION 

A century after its elaboration, the Drude formula is still 
commonly used to analyze the electrodynamic response of 
metals. 1 ' 2 However, the conceptual basis underlying Drude's 
phenomenological description has been entirely changed 
by the modern quantum theory of the solid-state. The key 
point is the emergence of long-lived quasiparticle excita- 
tions at low energy and low temperature, which are the 
actual charge carriers in good metals. Understanding the 
emergence of quasiparticles from a correlated liquid of in- 
teracting electrons is the great achievement of Landau's 
Fermi liquid (FL) theory, 3 which is precisely half as old as 
the Drude theory. In FL theory, the existence of long-lived 
quasiparticles is due to the vanishing of their scattering rate 
as the Fermi surface is approached, because of phase-space 
constraints and the Pauli principle. For quasiparticles to 
be well-defined, however, they must have a relaxation rate 
smaller than their energy (and hence than the available ther- 
mal excitation energy ~ fc B T), a condition which is met at 
low temperature close to the Fermi surface of simple metals. 

When exciting carriers at a low frequency co, one induces 
intra-band transitions of energy hco between states within 
windows of order k B T, thus probing the relaxation rate up 
to energies typically hco + k B T. A Drude-like response is ex- 
pected if the relaxation rate does not vary appreciably over 
this energy range. For Landau quasiparticles in a FL, the re- 
laxation rate increases quadratically with increasing energy 
and temperature. The conditions for a Drude response are 
therefore met if hco <C k B T, but deviations are expected if 
hco >k B T. Deviations from a pure Drude behavior are actu- 



ally common in metals, and this is often somewhat loosely 
interpreted as a violation of Fermi-liquid behavior. How- 
ever, in order to distinguish non-Drude features which are 
consistent with — or even signatures of — FL behavior, from 
those which indicate a genuine breakdown of FL theory, one 
needs to understand the response of Landau quasiparticles 
to electromagnetic waves. This is a difficult task in gen- 
eral, which has regained interest recently, 4 ' 5 along with the 
development of low-frequency spectrometers. 6,7 

In this article, we address this problem in the context of 
local Fermi liquids. By 'local', we mean that the scattering 
rate, and more generally the single-particle self-energy, is 
independent of momentum. From a theory viewpoint, local 
FLs are realized e.g. in the Dynamical Mean Field Theory 
(DMFT) of strongly correlated electron models. The success 
of this approach in practice also demonstrates that a weak 
momentum-dependence is a reasonable approximation for 
a wide class of correlated electron materials, at least in 
some range of composition, temperature, etc., when spatial 
correlations are short-ranged. 

The considerable simplification resulting from locality al- 
lows us to derive here analytically a universal scaling form 
of the optical conductivity in the FL regime. We show that 
a Drude-like behavior is indeed recovered in the regime 
hco <C k B T. In contrast, clear departure from Drude be- 
havior is found at higher frequency hco ~ k B T when the 
frequency dependence of the scattering rate becomes impor- 
tant. Characteristic non-Drude signatures of FL behavior in 
the optical conductivity are identified in this regime. 

We also perform DMFT calculations of the optical conduc- 
tivity for a microscopic model of a hole-doped Mott insulator. 
We show that the results accurately obey the FL scaling ex- 
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pressions, and that characteristic non-Drude features of the 
DMFT optical conductivity in the thermal regime hco ~k B T 
are explained by the FL scaling analysis. The model calcula- 
tion also allows to clearly identify the limitations and range 
of applicability of universal FL behavior. 

Finally, we discuss the conditions for a possible experi- 
mental observation of the FL universal scaling laws and FL 
signatures in optical measurements. 

The paper is organized as follows. Section II reviews 
the Drude theory and its generalizations. In Sec. Ill, we 
give the general formula for the optical conductivity of a 
system with a local self-energy. We then introduce the low- 
energy expression of the self-energy in a local FL, derive 
analytically a universal scaling expression for the optical 
conductivity, and discuss the different regimes of physical 
relevance. Section IV presents a comparison of the FL scaling 
laws with DMFT calculations. In Sec. V, we address issues 
related to the experimental observation of FL scaling laws 
in the optical conductivity, and discuss in more details the 
case of UPd 2 Al 3 . Our conclusions are given in Sec. VI, and a 
series of appendices collect additional material. 



II. PRELIMINARIES - DRUDE THEORY AND BEYOND 

In the Drude theory of conduction in metals, a single 
frequency-independent time t d governs the relaxation of 
the current. The assumption is that the current decays 
exponentially once the electric field has been turned off. The 
classical equation of motion then leads to the dc conductivity 
cr dc = ne 2 T D /m, with n the carrier density and m the carrier 
mass. In the presence of an oscillating electric field, the 
complex frequency-dependent conductivity cr(co) = a x (co) + 
icr 2 (a>) reads: 



o" dc 



m — ia)+l/T D 1 — ia>T D 



(1) 



In many materials however, especially those with strong 
electron correlations, a single frequency-independent re- 
laxation time does not provide an accurate description 
of the experimental data. 8 In order to describe the full 
frequency-dependence of the conductivity, it is convenient 
to introduce 9 a memory function M{co) such that: 



co + M(co) 



(2) 



In this expression, we define the plasma frequency from the 
sum-rule over the whole frequency range: 



dco cjjCco). 



(3) 



The complex function M(a>) = M^co) + iM 2 (co) plays the 
role of a self-energy for the optical response. The defini- 
tion (3) of co 2 p ensures that co dominates over M(co) in 



the expression of a{co) at large frequencies. Indeed, the 
Kramers-Kronig relations: 
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(4) 



imply that <J 2 (co ~ oo) ~ 2 f dQ cr 1 (r2)/(7ia)). 

Expression (2) can be cast in a form which is formally 
analogous to the Drude expression ("generalized Drude 
model"): 



a(co)- 



p _ 



ico + 1/t t (w) 



(5) 



with: 



1 + 



1 opt 



0) 



®(6>)M 2 (&>). (6) 



Note that the complex conductivity obeys c*(co) = <r(— co) 
under complex conjugation, so that cr 1 , M 2 , S> (resp. cr 2 , 
M x ) are even (resp. odd) functions of frequency. @(a>) is 
often denoted m/m*(co), hence defining an optical effec- 
tive mass, and T opt (co) is often written as t*(o>). At low 

frequency, @(w — » 0) = (l + S^Mj^o) -1 renormalizes 
the bare plasma frequency and the spectral weight of the 
Drude peak. In many cases, the frequency-dependence of 
the optical scattering rate l/T opt (a>) at low frequency is 
mainly determined by M 2 (co), with 3>(co) having a milder 
frequency dependence (see below). 

When analyzing experimental data, the imaginary part of 
the memory function which controls the optical scattering 
rate can for example be determined through: 



M 2 (co) = Re ■ 



(7) 



cf(co)' 

with the plasma frequency determined by the sum rule (3). 

in. OPTICAL CONDUCTIVITY OF LOCAL FERMI LIQUIDS 

We now specialize the discussion to the optical conduc- 
tivity of local Fermi liquids, i.e. systems in which the single- 
particle self-energy S obeys the low-frequency low temper- 
ature behavior of Landau Fermi-liquid theory and, addition- 
ally, does not depend on momentum. This is the case, in 
particular, of strongly correlated electron models and ma- 
terials treated in the framework of dynamical mean-field 
theory. 10 

We first recall the simplifications encountered in the Kubo 
formalism in this context, and the resulting expression of 
the optical conductivity. We then show that universal scaling 
laws emerge in the Fermi liquid regime and discuss these 
laws in the different physical regimes. 
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A. General expression of the optical conductivity 

When the single-particle self-energy E(e) has no momen- 
tum dependence, the Kubo formalism leads to the following 
general expression of the optical conductivity: 



d£[f[E)-f(£ + hco)] 
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d^^)A^,e)A^,e + hcj) (8a) 



in which /(e) is the Fermi function, while $(£;) and P$£, e) 
are respectively the transport function and the one-particle 
spectral function, defined in detail below. The derivation 
of this formula is outlined in Appendix A. The key point is 
that vertex corrections associated with the current-current 
correlation function exactly vanish in the case of a local 
(momentum independent) theory, 11 so that the conductivity 
can be entirely expressed in terms of the one-particle self- 
energy. 

Expression (8a) of the optical conductivity applies to a 
single-band system, to which our discussion is limited in this 
article for simplicity The entire information about the band 
dispersion is encoded in the transport function <J>, defined 
by: 



2e 2 



(8b) 



In this expression, L is the system size, d is the dimensional- 
ity) ?fc = E k — fx is the dispersion of the non-interacting 
Bloch band measured from the chemical potential, and 
v k = (l/ft)V fc £ fc is the corresponding group velocity. Many- 
body effects enter through the single-electron spectral func- 
tion, A(E,, e) which is related to the self-energy E(e) = 
S 1 (e) + iE 2 (e)by 



-E 2 (e)/7t 



2 ' 



(8c) 



In all these expressions, e /ft designates a frequency, while 
the momentum dependence of the spectral function stems 
from E, k . 

The transport function $(<!;) is usually a slow function 
of its argument, in contrast to A(E,, e) which in a Fermi 
liquid peaks at e E,. When the energy dependence of 
<£(<?) is negligible over the energy range where the spectral 
functions are appreciable, the second integral in Eq. (8a) 
reduces to the convolution of two Lorentzian functions, and 
the real and imaginary parts of the complex conductivity can 
be recast into the very convenient form 12 " 14 (Appendix A) : 



o-(w): 



;<Ko) 



de 



f{e)-f(e + hco) 
ha> + E*(e) - E(e + hco) ' 



(9) 



Equation (9) reduces to the Drude formula Eq. (1) if we 
put E(e) = — ift/(2T D ), and if 3>(0) is evaluated using the 
three-dimensional isotropic electron gas formulae. 



In the context of a single-band system, the total sum rule 
as defined from Eq. (3) reads: 



r °° 5$ 2e 2 



(n(£ fc )V 2 £ fc ) BZ (10) 



in which (• • • )bz designates an average over the Brillouin 
zone, and n(£ fc ) = (c k c k ) = J de f{e)A{k, e) is the momen- 
tum distribution of the electrons. For a simple tight binding 
band E k oc ^ a cosfc a , the right-hand side becomes propor- 
tional to the kinetic energy. Note that the value of the sum- 
rule thus depends on temperature and interaction strength. 
This is in contrast to the case where the whole solid (with 
all bands) is considered, in which case the plasma frequency 
co 2 is set by the bare electron mass and total number of 
electrons (/-sum rule). Note also that the approximate for- 
mula (9) in which the transport function was replaced by 
a constant does not reproduce correctly the total sum-rule, 
which is expected since it is only valid at low-energy. We 
shall see however that it can in general be used to reliably 
estimate the spectral weight of the Drude peak. 



B. Scattering in a local Fermi liquid 

The self-energy of a local Fermi liquid reads, at low fre- 
quency and temperature: 



E(e, T) = ( 1 - - | f - 



Znk B T 



[e 2 + (7ifc B T) 2 ] . (11) 



In this expression, Z is the quasiparticle spectral weight. In a 
local Fermi liquid, it is also related to the quasiparticle mass 
renormalization through Z = m/m* (i.e. the quasiparticle 
Fermi velocity is v* = Zv F ). Close to the Fermi surface 
and for e <C k B T, the spectral function is approximately 
a Lorentzian peak of weight Z, centered around e = ZE, k 
(the quasiparticle dispersion) . We define the temperature- 
dependent quasiparticle lifetime from the zero-frequency 
value of E 2 as: 



ft (fc B T) 2 
— = 2Z|E 2 0? = 0, r)| = 2ti-^— 



(12) 



-qp 



It corresponds to twice the width of the Lorentzian quasi- 
particle peak in the spectral function, that is to the decay 
rate of the probability (square of the Green's function). 

Fermi liquid behavior self-consistently relies on the exis- 
tence of a sharp Fermi surface, and applies when the scat- 
tering rate of the typical excitation is smaller than its energy. 
The scattering rate ft/T qp must be compared to a typical 
thermal excitation energy, say 2nk B T. The characteristic 
Fermi liquid energy scale fc B T is defined as the temperature 
where fi/i^ = 2nk B T, and will serve as a basic unit below. 
This is illustrated in Fig. 1 (left panel) which sketches the 
dependence of the quasiparticle lifetime on temperature. 
The temperature T Q also sets the coupling strength: a small 
T corresponds to strong electron-electron interactions and 
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FIG. 1. Fermi liquid scattering rate ftf^ 1 . (a) At zero energy, Trr^ 1 
increases quadratically with temperature (solid line) . The dashed 
line indicates 2nk B T. The intercept defines a temperature scale 
T and a frequency scale hco = 2nk B T . (b) At finite T, ftT _1 (w) 
increases quadratically with the frequency co. The applicability of 
the frequency-independent (Drude) approximation is limited to 
low frequencies. 



a large value of E 2 . Hence, the quasiparticle lifetime dimin- 
ishes rapidly as a function of temperature and frequency. 
The precise relation of the scale T to the scale T FL , be- 
low which strict Fermi liquid behavior holds, is discussed 
in Sec. IV. In strongly correlated (local) Fermi liquids in 
which Z is small, e.g. close to a Mott transition (Sec. IV) 
or in heavy-fermion materials, all physical quantities scale 
with a single energy scale, and T is proportional to ~ ZD 
(with possibly a small value of the prefactor) . Here, D is the 
half-bandwidth of the bare band. The factor 1 /Z which has 
been pulled out in front of E 2 in order to define T ensures 
that in such a case E 2 is a scaling function of T/ZD and 
e/ ZD. 

The right panel of Fig. 1 suggests that the frequency- 
dependence of the quasiparticle scattering rate is actually 
important. A Drude-like optical response with a constant re- 
laxation time t d ~ Tqp can only be expected to provide a rea- 
sonable description in the very low-frequency or relatively 
high-temperature regime hco < 2nk B T. Whenfto) > 2nk B T, 
the energy dependence of S(e) cannot be neglected and a 
non-Drude response arises, as discussed in the following 
section. 

We finally note that a real, frequency-independent 
Hartree term E^O, T) should in fact be added to Eq. (11) 
for completeness. It sets the location of the Fermi surface 
from E k = [I — E^O, T) and can be viewed as a shift of 
the chemical potential, which will be omitted for simplic- 
ity in all equations. An imaginary frequency-independent 
part can also be added, to mimic the effects of the impurity 
scattering. This is considered in Appendix E. 



C. Scaling form of the optical conductivity in a local Fermi 
liquid 

We now show that the optical conductivity obeys a uni- 
versal scaling form in terms of the two variables oor^ and 
hco/(2nk B T). Whereas it reduces essentially to the Drude 
form in the low-frequency limit, its full frequency depen- 



dence is markedly different. 

The universal scaling form is derived by inserting Eq. (11) 
into Eq. (9). The calculations can be performed analytically 
and yield: 



o-(<o) = o- dc 5" 

.tt2 



hco 



2nk B T 



, COT 



qp 



(13a) 



'dc 



S?{x,y) 
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Z$(0)t 



qp 



du- 



(13b) 

1 + x 2 — iy + u 2 



The scaling function 5/ is evaluated and displayed in Ap- 
pendix B. One obtains: 



5"(x,y) 



6i 



n 2 xr{x,y) 



fyQ [l + r(x,y)-ix]) 



-V>Q [l + r(x,y) + ix])}, (13c) 

where r(x, y) = \J\ + x 2 — iy and ip is the digamma func- 
tion defined as i/>(z) = lim M _ 00 [logM - Y,n=o l/(" + z )] • 
Equation (13a) emphasizes the emergence of two natural 
time/ frequency scales: the quasiparticle time T qp , and a "co- 
herence" time h/(2nk B T). Alternatively, one can reexpress 



cr(w) = cr dc 5" 



T T< 



(14) 



with co = hco/(2nk B T ) and T = T/T Q dimensionless vari- 
ables normalized to the basic scale T . This emphasizes that 
the optical conductivity in a Fermi liquid is a scaling function 
of co/T and co/T 2 . Equations (13) constitute the key ana- 
lytical result of this article. They replace the Drude formula 
by a universal scaling form of these two frequency scales, 
which is valid for local Fermi liquids. Let us emphasize that 
all high-energy scales such as the bare bandwidth or the 
plasma frequency co p have disappeared from the scaling ex- 
pression (13a). Instead, only low-energy scales appear, such 
as T and Z$(0) (the latter is shown below to be related to 
the low-energy Drude weight) . 

In Fig. 2, we plot the real and the imaginary parts of 
<j(co)/<j dc as a function of hco/(2nk B T) on a log-log scale 
at a given low temperature T/T = 0.1. (The full fre- 
quency and temperature dependence is shown on a three- 
dimensional plot in Appendix B.) Three frequency regimes 
can be identified from this plot: 

• At low-frequency co < co L < 2nk B T /h (with co L to be 
made precise below, of order rZ} at low temperature) 
the conductivity follows closely the Drude model, with 
a saturation of a 1 below the characteristic frequency 
t^ 1 and a 1/co 2 decay above it. t^ 1 also separates the 
dissipative regime (with larger crj) from the inductive 
regime (with larger <t 2 ). 

• When co approaches 2nk B T/h, the conductivity de- 
viates from the Drude behavior: tr l displays a 
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10 u 

ftto/(2nk B T) 

FIG. 2. Fermi-liquid conductivity Eqs. (13) at low temperature 
(solid lines), below the temperature defined in Fig. 3. The blue 
(red) lines show the real (imaginary) part of the conductivity. The 
dashed lines show the low-frequency Drude-like behavior given 
by Eq. (15). The dotted lines show Eq. (18). The characteristic 
frequency scales co L and co H correspond to those defined in Fig. 3 
and Eq. (20). 



pronounced shoulder with much weaker frequency 
dependence — the feature appears as a shoulder in a 
log-log plot, as a 'foot' in a lin-lin plot, see below. In 
this 'thermal' regime, the conductivity behaves induc- 
tively (cr 2 > cjj) rather than dissipatively. 

• Increasing co further, a 2 decays more rapidly, and at 
o> H (to be defined below) the data become dissipative- 
like again (<j 1 /<j 2 > 1). This is actually a consequence 
of the assumed unbounded quadratic increase of the 
scattering rate with frequency, and might not be physi- 
cal in this already high frequency regime (see Sec. IV). 
Above co H , a 1 recovers a Drude-like l/co 2 decay, while 
a 2 turns to l/co 3 . The l/co 3 behavior is an artifact of 
extending the co 2 in the self-energy Eq. (11) to high 
energies. 15 

These different regimes, as well as the precise conditions 
determining the characteristic frequencies co LH are repre- 
sented on Fig. 3. We now discuss these different regimes in 
more details. 



D. Drude regime 

The lowest frequency regime can legitimately be called 
'Drude regime', since there the frequency-dependence of the 
scattering rate can be ignored. In this regime, the frequency 
hco is small with respect to 2nk B T, the inverse quasiparticle 
time t^ 1 is smaller than, or comparable to the tempera- 
ture, but the product coT qp = y can take arbitrary values 
(Fig. 3). This regime is thus described by taking the limit 
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FIG. 3. Regimes of conduction. The Drude regime applies if 
co < 2nk B T/h, and the thermal regime if co > 2nk B T /ft. The 
conductivity along the two horizontal cuts at fixed temperature is 
displayed in Figs. 2 and 6. 



x = hco/(27ik B T) — » in the scaling form (13c). This yields: 



C7 dc 7T Z 



4^ 



(15) 



1WT 



qp 



where xp' is the derivative of the digamma function. This 
is shown as the dashed lines in Fig. 2. Interestingly, this 
expression differs from the simple Drude form a[co)/a dc = 
1/(1 — iwT D ). However, as detailed in Appendix C, the 
frequency dependence of the exact expression (15) can be 
rather accurately approximated by the simple Drude form, 
provided t d is chosen in an appropriate manner (which is 
such that t d differs from Tqp however) . 

The spectral weight in the Drude peak can be estimated 



as: 



2 

n 



dcoa^co) = Z$(0) 



(16) 



in which co D is a cutoff defining the Drude regime. This 
expression can be established in two ways. The first is to 
perform a direct frequency integration of the scaling expres- 
sion (15) over the whole frequency range (i.e. for the scaling 
variable y = coz^ varying from to oo). Alternatively, since 
<3>(if;) can be taken as constant in the low-frequency range, 
one can use expression (9) and observe that it yields a high- 
frequency behavior of u 2 ~ Z$(0)/co. By Kramers-Kronig, 
this leads to Eq. (16). 

Two remarks are in order regarding expression (16). First, 
the Drude weight scales with the quasiparticle weight Z. 
Close to a Mott transition, for example, Z is expected to 
vanish and so does the Drude weight. Second, it should be 
noted that Z measures the renormalization of the Drude 
weight as compared to the non-interacting (band) value 
$(0) and not the spectral weight of the Drude peak relative 
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to the total integrated spectral weight e to 2 . The latter 
relative weight is given by Z<f>(0)/e to 2 & @(0) with 
given by Eq. (10). 



E. Thermal regime and emergence of the 2n factor in the 
optical scattering rate 



In the thermal regime, hto is comparable to 2nk B T and 
wT qp ;> 1. This corresponds to the limit of fixed x = 
hto/(2nk B T) and large y = toi^ ;> x. The scaling function 
Sf{x,y) has the following expansion for y —> oo: 



Sf(x,y -» oo) : 



12i 

n 2 y 



+ 



16 

2 2 



(l + x 2 ). 



(17) 



Using this expression and performing a large-y expansion of 
1/(7 = l/[cj dc 5^(x,y)], one directly obtains a generalized 
Drude form for the optical conductivity: 



cr(to) i 



Z<t>(0) 



-fed 4- l/T opt (w) 
with the optical scattering rate: 
ft 2 



(18) 



t t (w) 37lfc B T 



[(ftw) 2 + (27rfc B T) 2 ] . (19) 



The same result is found by approximating in Eq. (5) 
by @(0) w Z$(0)/(e co 2 ), identifying with Eqs. (13), and 
expanding for large y. Eq. (18) is displayed in Fig. 2 as the 
dotted lines. While the Fermi-liquid form of the quasiparticle 
lifetime (a one-particle quantity) involves (fto>) 2 + (7i:fc B r) 2 , 
it should be emphasized that the optical scattering rate 
(a two-particle quantity) involves instead a factor 2n in 
the combination (hco) 2 + (2nk B T) 2 . This seems to have 
been first emphasized by Gurzhi. 16 Experimentally, the 2n 
factor has not been observed so far in simple metals. As 
discussed in details in Sec. V, in several correlated metals 
a scaling of the optical scattering rate following {hoi) 2 + 
(pnk B T) 2 was reported, with values of p ranging from 1 to 
2_4_i7-22 d e p ar t; Ure f r0 m the Fermi-liquid value p = 2 
may be attributed to scattering mechanisms with a different 
frequency/temperature dependence. 4,5 In Appendix E, we 
show that a frequency-independent scattering rate does not 
change the value p = 2. 

At low temperature, Eqs. (18) and (19) allow to define 
two characteristic frequencies to L H from the condition to = 
T~p t (w), or equivalently a^to) = a 2 (to): 



4 



7lfc B T 




(20) 



Plotted as a function of T, these two frequencies form a 
dome in the (to, T) plane, defining the thermal regime 
(green line in Fig. 3) . Below this dome, the conductivity be- 
haves inductively. While the scales to LH can be determined 




FIG. 4. Ac resistivity 1/ct 1 (cl), 7) as a function of T for several 
frequencies a>. These curves correspond to vertical slices in Fig. 3, 
and the green bars indicate the temperature at which the dome is 
crossed for each frequency. The temperatures are measured in units 
of T , and the frequencies in units of a> = 2nk B T /h. The normal- 
ization of the resistivity is <r dc (r/r ) 2 = 7iftZ$(0)/(24fc B r ). 



from the crossings of a 1 and a 2 plotted as a function of 
frequency (Fig. 2), the crossing of the dome is most easily 
identified by a minimum in a plot of 1/a^to, T) as a func- 
tion of T, at a fixed finite frequency. This is illustrated in 
Fig. 4. As this determination of to L H only relies on a l , it 
might be the most direct way of checking FL behavior in 
experimental datasets. 

Finally, we plot in Fig. 5 the imaginary part of the memory 
function, defined according to Eq. (2) with the conductivity 
given by the Fermi-liquid expressions (13). Fig. 5a displays 
the crossover from the Drude to the thermal regime. M 2 (to) 
increases quadratically for to — » 0, with a temperature- 
dependent curvature given in Appendix D. In the thermal 
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FIG. 5. Imaginary part of the memory function in a Fermi liquid at 
T = 0.02T (thick red lines), (a) At frequencies to < t" 1 , M 2 (co) 
increases as to 2 (thin blue line), with a temperature-dependent cur- 
vature given by Eq. (D2). (b) In the thermal regime, M 2 ((x>)/M 2 (0) 
scales as (7r/3) 2 (l + x 2 ) with x = HtL>/(2nk B T) (dashed line). 
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regime, the expansion (17) leads to the following form: 



f e n o) 2 
" p 

Z$(0) 



+ t 



Z$(0) 37rftfc B T c 



[(ft6j) 2 + (2 7 ifc B T) 2 ] . (21) 



The imaginary part is shown in Fig. 5b. At co = 0, Eq. (21) 
extrapolates to a value larger than the exact value iM 2 (0) 
by a factor (n/3) 2 (see Appendix D). Similarly, the approxi- 
mations (18) and (19) for the conductivity deviate slightly 
from the exact conductivity in the limit co — » 0. Setting 
co = in these equations yields cr(0) = (3/7i) 2 cr dc . All these 
observations are due to the fact that the expressions (19) 
and (21) holds in the thermal regime, but are not accurate 
at very low frequency. 



F. Coherent and "incoherent" regimes 
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FIG. 6. Fermi-liquid conductivity Eqs. (13) at high temperature 
(solid lines), above the temperature T 1 defined in Fig. 3. The blue 
(red) lines show the real (imaginary) part of the conductivity. The 
dashed lines show the low-frequency Drude-like behavior given by 
Eq. (15). The dotted lines show Eq. (18). 



It is seen from Fig. 3 that the three frequency scales co L H 
and 2nk B T/h defined above merge when the temperature is 
raised above T x = 3T /8. Above 7\, an "incoherent" regime 
is found in which the thermal scale 2nk B T /ft becomes the 
only characteristic frequency scale and cr 2 becomes smaller 
than a i at all frequencies, because t"* > co (see Fig. 6). 
The term incoherent is put in quotes, because at such a high 
temperature T ~ T the scaling form of the self-energy may 
be no longer valid. The actual incoherent regime occurs 
due to the breakdown of the Fermi liquid form Eq. (11) 
altogether. The question whether the scaling form applies is 
then irrelevant. 

The reciprocal argument, on the other hand, is valid. 
Observing two frequencies at which <j 1 = a 2 indicates that 
the Fermi liquid range has been reached. On a log-log plot 
of cocji^co^ versus co, these two crossing points coincide 
with two symmetric maxima located at co L H , separated by 
a minimum at co = 2nk B T /ft. 



IV OPTICAL CONDUCTIVITY OF A HOLE-DOPED MOTT 
INSULATOR WITHIN DYNAMICAL MEAN-FIELD THEORY 

A. Model and methods 

In this section, we present calculations of the optical con- 
ductivity for a specific microscopic model, the single band 
Hubbard model of a hole-doped Mott insulator. The calcu- 
lations are performed within single-site dynamical mean- 
field theory (DMFT). 23 As we shall see, the universal scaling 
form derived above allows one to identify specific signatures 
of Fermi-liquid behavior in the DMFT optical conductivity, 
which have not been emphasized previously. Conversely, the 
model calculation allows for a test of the scaling theory and 
especially of its range of validity as a function of frequency 
and temperature. Optical conductivity has been calculated 
with DMFT by several authors in various contexts (see e.g. 



Refs. 24 and 25, and, for a review and more references, 
Ref. 8). However, an explicit analysis in connexion with FL 
scaling laws has not been made, and requires calculations 
with very high accuracy solvers at low energy, which became 
available only recently. 

A semicircular density of states with half-bandwidth D has 
been used, and the model is considered in its paramagnetic 
phase. The DMFT equations were solved using highly accu- 
rate continuous-time quantum Monte Carlo (CT-QMC) 26,27 
and numerical renormalization group 28 ' 29 (NRG) solvers. 
The imaginary-frequency data from Monte Carlo were an- 
alytically continued using Pade approximants. Excellent 
quantitative agreement between the data obtained using 
the two techniques was achieved in a broad temperature 
and frequency range. At the lowest temperatures (below 
k B T/D = 0.01) and the lowest frequencies, however, arti- 
facts associated with the discretization of the energy mesh 
and spectral broadening become visible in the NRG data, 
and the CT-QMC+Pade method becomes preferable. Con- 
versely, the CT-QMC+Pade data becomes less accurate at 
higher temperatures k B T/D > 0.05, especially for larger 
frequencies hco/D > 1. The data used here are obtained 
by taking the low frequency part (hco/D < 0.15) from the 
calculation based on the CT-QMC+Pade self-energies, and 
the high-frequency part (hco/D > 0.15) from the NRG self- 
energies, with perfect matching in the intermediate region. 



B. Self-energy, local Fermi-liquid behavior, and key 
temperature scales 

The analysis presented here is based on the DMFT dataset 
used earlier in Ref. 30. There, transport and thermodynamic 
properties were discussed in detail, but only some aspects 
of the optical conductivity were addressed. For convenience, 
we restate here the key temperature scales identified in this 
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an order of magnitude smaller. Therefore, when analyzing 
the DMFT results in the light of Fig. 3 and of the scaling 
analysis, it should be remembered that FL behavior actually 
fully applies only below T FL = 0.1T . 

At higher temperatures, deviations from FL rapidly appear 
for electron-like (co > 0) single-particle excitations (Fig. 7). 
There, the deviations from the parabolic form become sub- 
stantial at a frequency hco + 7iT FL , where the real part 
of the self-energy (not shown, see Ref. 30) displays a kink. 
For hole-like excitations (co < 0), the parabolic behavior 
is more robust, and the kink appears only at co_ as 0.2D. 
The transport (resistivity, thermopower) probes a frequency 
window of a few k B T, and these quantities deviate from the 
FL universal behavior when k B T > hco + . 



FIG. 7. DMFT self-energy scaling plot. At low temperatures, the 
curves collapse to a parabola. By comparing with Eq. (11) and 
taking into account that one has Z = 0.22 for a doping 5 = 0.2, 
one can determine k B T « 0.1D p» 10fc B T FL . 



previous work, as well as their evolution with doping. The 
actual data considered in detail later in the present paper 
are for a doping level 5 = 20% and a coupling U/D = 4 
(at which the undoped system is a Mott insulator), but we 
mention the doping evolution of the key quantities. 

The quasiparticle weight Z was found to be approximately 
equal to the doping level Z f*s 5 (more precisely, Z = 0.22 
was found for 5 = 0.2). A temperature scale of key impor- 
tance is the Fermi-liquid temperature scale T FL . This was 
defined as the temperature below which the T 2 behavior 
of the resistivity and the co/T scaling of the self-energy 
apply. From the DMFT data, this scale was identified as 
TpL = 0.055D, which gives 0.01D for 5 = 0.2, a very low en- 
ergy scale. In the same work, 30 the quasiparticle features in 
the spectra were shown to persist to a much higher temper- 
ature, 5D. Important particle-hole asymmetry was found in 
many physical properties. At a scale approximately equal to 
5D, these 'resilient quasiparticles' disappear and the system 
becomes an incoherent 'bad metal' with a resistivity which 
becomes larger than the Mott-Ioffe-Regel value. It should 
be emphasized that the Fermi-liquid scale T FL is one order 
of magnitude smaller than the scale at which the crossover 
into the bad-metal regime occurs. The latter corresponds 
quite accurately to the Brinkman-Rice scale ~ 5D, which 
corresponds to the renormalized kinetic energy of the quasi- 
particles and is much larger than T FL (although both scales 
are proportional to the doping level). 

In Fig. 7, we display the imaginary part of the self- 
energies S, as a scaling plot -ImS(w, T)D/(nk B T) 2 versus 
hco/(nk B T). The data nicely obey the Fermi-liquid law 
oc 1 + (hco/nk B T) 2 at low temperatures. By comparing 
with the prefactor of this scaling behavior as defined in 
Eq. (11), one determines the scale T to be (at 5 = 20%): 
k B T » 0.1D, so that T » 10T FL . For an arbitrary doping 
level, one finds, k B T w 0.575D. The scale T Q is thus rather 
close in magnitude to the Brinkman-Rice scale, while T FL is 



C. Optical conductivity at low temperature 

Fig. 8 for T = T FL /4 = T /40 demonstrates that the 
optical conductivity of the hole-doped Hubbard model ob- 
tained from DMFT is very well described by the universal 
FL scaling form (13) derived in the previous section, in the 
low-frequency and low-temperature regime. This is indeed 
expected from the previous figure demonstrating FL scaling 
of the single-particle self-energy. 

When looking at cr^co) on a lin-lin scale (Fig. 8a), the 
narrow Drude peak at low-frequency is followed at higher 
frequency hco ~ 2nk B T by a characteristic non-Drude 'foot'. 
On a log-log plot (Fig. 8b), this appears as a shoulder. This 
does not signal non-Fermi liquid physics, as one might 
naively think, but is actually a key feature of FL theory which 
signals the onset of the thermal regime at hco ~ 2jzk B T. 
The non-Drude foot/ shoulder is thus, somewhat counterin- 
tuitively, a striking signature of Fermi-liquid behavior. 

At higher frequencies, the DMFT data displays two peaks. 
The first one, for hco w 0.5D (corresponding typically to the 




Hoj/D hco ID 

FIG. 8. Comparison of the optical conductivity at k B T/D = 
0.0025D (T = r FL /4) for the doped Hubbard model calculated 
within DMFT (solid lines) to the universal FL scaling form (dashed 
lines). The real (blue) and imaginary (red) parts are plotted on a 
lin-lin (a) and log-log (b) scale. 
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mid-infrared MIR regime) is associated with the transitions 
between the quasiparticle band and the lower Hubbard 
band. The high-frequency peak at hco fa AD ~ U corre- 
sponds to the transitions to the upper Hubbard band. These 
peaks are, obviously, not present in the FL expressions. Like- 
wise, at the highest frequencies the bare particle dynamics 
(with <j 2 (X is recovered in the DMFT data, whereas 
extrapolating FL behavior to infinite frequencies would lead 
to the incorrect behavior cr 2 oc l/<w 3 . 



D. Memory function 

More subtle corrections to the FL are seen in the DMFT 
data, also at frequencies smaller than the MIR, in the 
'foot/ shoulder' region (but above the low-frequency kink of 
the self-energy) . To resolve them more clearly, it is conve- 
nient to look at the memory function. In Fig. 9 we plot the 
memory function at a low temperature T = T FL /4. 

The memory function has features on crossing the Hub- 
bard bands (Fig. 9a), which we will not discuss here. Below 
the MIR scale (Fig. 9b), approximately linear and quadratic 
behaviors are seen for the real and imaginary parts, respec- 
tively (Fig. 9c). On zooming up further (Fig. 9d), one sees 
an excellent agreement between the DMFT results and the 
FL scaling form, including the small crossover at the lowest 
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FIG. 9. Memory function at k B T/D = 0.0025 (T = T FL /4). (a)- 
(d): real (solid blue) and imaginary (solid red) parts of the memory 
function, compared with the FL scaling forms (dashed) . The data 
are shown for several frequency windows ranging from a very 
broad one (a), to a very narrow one (d). Panels (e) and (f) show 
the frequency derivative of M 1 (solid lines) compared with the FL 
expressions (dashed lines) . 



frequencies co = I/t™ 0.001. 

Whereas agreement between the DMFT data and the FL 
scaling form is perfect at the lowest frequencies, some devi- 
ations appear at a small but well-defined and important FL 
frequency scale, the frequency to + 0.03 associated with 
the positive-energy (electron-like) 'kink'. This is seen clearly 
in M 2 or the derivative dM 1 /dco (Fig. 9e and f)- The devia- 
tion from the FL form goes into the direction of a smaller 
memory function. This is due to the fact that, in the hole- 
doped Hubbard model studied here, the positive frequency 
'resilient quasiparticles' scatter less than the parabolic be- 
havior from Landau FL theory would predict. We stress that, 
for the hole-doped Hubbard model, the corresponding de- 
viation cannot be observed in photoemission spectroscopy, 
since they appear at positive excitation energies. Optical 
spectroscopy, because it probes particle-hole excitations, is 
thus a powerful tool which could probe the existence of 
these resilient quasiparticles. Moreover, combining photoe- 
mission with precise optical spectroscopy could be used to 
reveal the strong particle-hole asymmetry of these excita- 
tions. Signatures of the low-frequency kink from optical 
spectroscopy would also be very interesting. 



E. Temperature-dependence of the optical conductivity 

Finally we discuss the evolution of the optical conduc- 
tivity at higher temperatures. Fig. 10 displays cr 1 (<x)) for 
several temperatures (D/k B T = 400,200,100,50,20,10). 
The DMFT data are compared to the analytical form dic- 
tated by the FL scaling function, in which is determined 
from the FL T 2 dependence. The difference between the 
FL and the DMFT becomes pronounced for T > T FL . Note 
that a milder frequency dependence is seen above T FL , both 
in the FL and DMFT (in the former it occurs due to the 
narrowing of the thermal crossover foot/shoulder). This 
warns again that interpreting apparent non-Drude power 
laws as a signature of non-Fermi liquid behavior is a risky 
enterprise. 



10° pr 



10" 



10" 



10" 



10" 



10- 




— k B 770 = 0.1 

fc B r/D = 0.05 

k B T/D = 0.02 

— fc B r/D = o.oi 

k B T/D = 0.005 

k B T/D = 0.0025 

mi I 



10 



io- 2 io- 1 

hco/D 



FIG. 10. Comparison of the real part of the optical conductivity 
in the Hubbard model from DMFT (solid lines) with the analytical 
FL scaling expressions (dashed lines). Inverse temperatures are 
D/k B T = 400, 200, 100, 50, 20, 10. 
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FIG. 11. (a) Temperature dependence of the imaginary part of 
the memory function in DMFT. The temperatures and the color 
code are the same as in Fig. 10. (b) Scaling plot. The ha>/(2nk B T) 
scaling applies at low temperatures and intermediate frequencies. 
The dashed line shows the parabolic behavior corresponding to 
Eq. (21). 



In Fig. 11, we show the temperature evolution of the 
imaginary part of the memory function. In Fig. lib the 
data are plotted as M 2 /T 2 versus co/T. This reveals clearly 
the scaling behavior, consistent with the (ftw) 2 + (27rfc B T) 2 
dependence in the thermal regime of FL theory. 

At low T, FL deviations from the quadratic dependence 
on frequency are seen at the lowest frequencies, indicating 
the onset of the Drude regime where hco <C 2jzk B T. On 
the high energy side, the discrepancies are seen above oj + . 
Above TpL the data deviates from a parabola and the co/T 
FL scaling does not apply anymore. 



V IMPLICATIONS FOR EXPERIMENTS 

In this section, we address the implications for the exper- 
imental optical signatures of Fermi-liquid behavior, in the 
light of the theoretical results discussed in Sec. III. We will 
also discuss what has been obtained until now in a num- 
ber of materials, in particular the heavy-fermion materials 
CePd 3 , UPd 2 Al 3 , and URu 2 Si 2 , heli-magnetic MnSi, the or- 
ganic conductor K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br ;c Cl 1 _ x , the 
doped semiconductor SrTi0 3 , and the doped Mott-Hubbard 
insulator HgBa 2 Cu0 4+5 . We will see that the experimental 
optical data published until now do not provide a sufficiently 
broad spectral range to distinguish all of the optical features 
of a Fermi liquid. The first signature is the 'foot/shoulder', 
illustrated in Figs. 2 and 8, which marks the deviation from 
a low-frequency Drude-like behavior. This feature occurs at 
the frequency co ~ 2nk B T /h — in-between two frequencies 
where the imaginary and real parts of the conductivity are 
equal — and disappears as the temperature is raised above 
T 1 (Figs. 3 and 6) . The second signature is the characteristic 
frequency/temperature scaling of the optical relaxation rate 
shown in Eq. (19), or the equivalent behavior of the memory 
function given in Eq. (21). While the first signature can in 
principle be observed in the raw data, the second requires a 
determination of the spectral weight, in order to invert the 
complex conductivity. A third signature — present in the raw 
data — is the minimum in l/cr 1 , when plotted as a function 
of temperature at finite frequency (Fig. 4). 

Experimentally, the most direct clue of Fermi-liquid 



physics is a T 2 law in the resistivity. This is usually observed 
between a ground-state ordering temperature T c , and a 
scale T FL above which additional scattering mechanisms 
contribute to the resistivity with a different temperature 
dependence. For T c < T < T FL , the Fermi-liquid signatures 
are expected in the conductivity (provided T < 7\ = 3T /8) 
around the frequency v = k B T /ft, which is ~ 1.3 THz for 
T = 10 K, corresponding to 44 cm -1 . The signature may be 
masked by optical phonons — with energies typically above 
5 mey i.e. 1 THz or 40 cm -1 — and inter-band transitions. 
Another possible limitation is impurity scattering, which 
reduces the value of T 1 as discussed in Appendix E. 

Fermi-liquid behavior of the dc resistivity has been re- 
ported in a variety of materials. Using a bolometric direct 
absorption technique, Webb et al. 31 observed the narrow 
zero-frequency mode corresponding to the Drude peak for 
CePd 3 . The presence of the peak was also indicated by 
Awashti et al. 32 using resonant cavities for three discrete 
frequencies (for an extensive review of optical properties of 
heavy-electron compounds, including Fermi-liquid aspects, 
see Ref. 33). 

UPd 2 Al 3 is a heavy-fermion material with a mass enhance- 
ment of 66. 34 The 5/ moments order anti-ferromagnetically 
below 14.3 K, 35 and a superconducting phase coexisting 
with the magnetic order develops below 2 K. In the temper- 
ature range between these two transitions, the resistivity in- 
creases as T 2 with a typical coefficient A ~ 1 /ifi cm/K 2 , sug- 
gesting Fermi-liquid behavior of the charge carriers 36,37 with 
a TpL of the order of 15 K. For this material, o{oS) has been 
measured in the range from 0.002 to 1.3 cm -1 using a coax- 
ial technique in Corbino geometry. 38 Figure 12 shows the mi- 
crowave data of Ref. 38, measured at T = 2.75 K, well below 
T FL . The conductivity can be well fitted to the Drude model 
Eq. (1), with the parameters a dc = 0.105 (/xfi cm) -1 and 
t d = 4.7 x 10~ n s. This led the authors to the conclusion 
that impurity scattering dominates in this frequency range. 




Frequency v (GHz) 

FIG. 12. Microwave optical response of UPd 2 Al 3 at T = 2.75 K. 
The thin lines show the measurements of Ref. 38. The solid lines 
show a fit to Eqs. (13) with the parameters T = 350 K and 
Z$(0) = e (1.6x 10 14 s -1 ) 2 . Inset: same data on a wider frequency 
range. The dashed lines show the Drude model. 
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However, the extrapolated dc resistivity of ~ 9 /xfi cm is 
of the order AT 2 , which may also indicate that electron- 
electron scattering dominates the dc resistivity. As a matter 
of fact, the data can be equally well fitted to the Fermi-liquid 
formula (13), as show in Fig. 12. This provides an alterna- 
tive interpretation for the success of the Drude model: the 
scattering rate is due to electron-electron interactions, but 
it is dominated by the temperature in this low-frequency 
range. At the experimental temperature T = 2.75 K, the 
thermal regime where deviations from the Drude model 
due to the frequency dependence of the scattering rate are 
expected, is around k B T/h = 360 GHz (inset of Fig. 12). 
Note that again T FL is an order of magnitude smaller than 
the value T Q ~ 350 K resulting from the fit in Fig. 12. 

Recently, Nagel et al. 22 reported optical data for URu 2 Si 2 , 
exhibiting the l/T opt (w, T) oc (hco) 2 + (pnk B T) 2 depen- 
dence of the optical scattering rate with p & 1, i.e. 
well below p = 2 expected for a Fermi-liquid (Eq. 19). 
To our knowledge, the two cases in which a value 
of p closest to 2 has been observed are the zc-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br x Cl 1 _ x organic compound (p as 2.4) 39 
and the underdoped HgBa 2 Cu0 4+5 (p & 1.5). 21 The optical 
conductivity of K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br x Cl 1 _ x does 
have a narrow Drude peak followed by a foot. However, the 
frequency region co < co L has not been fully explored. 

For most Fourier-transform spectrometers, the lower limit 
of the spectrometer range is about 20-40 cm -1 . As a result, 
experiments typically see only the upper frequency part of 
the shape shown in Fig. 2, containing the characteristic fre- 
quency scales co H , while the part of the spectrum around co L 
is usually not reported. This upper frequency part is itself 
often approximated by a Drude form, which however should 
not be confused with the true low-frequency behavior of the 
conductivity. For a number of materials, the experimental 
reports have hinted toward a cross-over to a low frequency 
regime which could be Fermi-liquid like. For example, in 
MnSi the resistivity is found to follow p = AT 2 in the heli- 
magnetic phase below 30 K, while the optical conductivity 
is given by ct^co) oc o) -0 5 down to 30 cm -1 , even for data 
taken at 10 K, indicating that a cross-over to Fermi-liquid 
like behavior, if present, would have to occur below 30 cm -1 
(Ref. 40). 

This is remedied to a large extent by state of the art time- 
domain THz spectrometers, which typically span a frequency 
range from 1 up to 100 cm -1 . The limiting factor in the 
latter case is provided by diffraction in the long-wavelength 
limit, imposed by the finite sample size which is typically of 
the order of a few mm 2 . In Nb-doped SrTi0 3 , T 2 resistivity 
in a broad temperature range is an indication for Fermi- 
liquid behavior. 41,42 Time-domain THz data reveal a very 
narrow zero-frequency mode, 43 similar to that in Fig. 2. On 
the other hand, whether or not the foot/ shoulder at higher 
frequency is present in the data, has been impossible to 
establish because strong optical phonons mask the electronic 
part of the spectrum in the relevant frequency range. 

A team involving the present authors recently reported 
Fermi-liquid like features for underdoped high-T c cuprates 
in the pseudo-gap phase, in particular the scaling collapse 



of the form l/z opt (co, T) oc {hoi) 2 + (pnk B T) 2 for a broad 
range of co and T. 21 Further analysis of these data in relation 
to the features discussed in the previous sections is under 
way. 

The examples discussed in this section illustrate the exper- 
imental challenges that need to be met in order to establish 
in a single material the different aspects of the Fermi-liquid 
optical conductivity, such as displayed in Fig. 2. To the best 
of our knowledge, this is not available in the literature. Still, 
with appropriate choice of material, sample geometry and 
optical instrumentation, those challenges may ultimately be 
met. 



VI. CONCLUSION 

We have derived analytical formula for the universal scal- 
ing laws which describe the optical response of local Fermi 
liquids. These laws depend on two variables, ftw/(27ifc B T) 
and wTqp = hco/(2nk B TXT Q /T), and reveal clear signa- 
tures of FL behavior, which are not contained in the simple 
Drude form. These new laws can help distinguishing, among 
the ubiquitous deviations from Drude-like behavior, those 
which prove FL behavior from those which disprove it. We 
have illustrated this by comparing the scaling forms with 
DMFT data for the Hubbard model at low doping. Finally, 
we have reviewed a number of experimental works, and con- 
cluded that several of these characteristic features remain 
to be observed. 
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Appendix A: Kubo formula for a local Fermi liquid 

In general, the long-wavelength linear conductivity tensor 
is related to the current susceptibilities by the Kubo formula 



CT a/3 (w) = — [Xafs(<0) - XapW] ■ (Al) 



X a p{oj) is the retarded macroscopic current-current corre- 
lation function, which may be obtained from the analytic 
continuation of the corresponding imaginary-time function: 



Xa[t(\n n ) = - — 



dze in ^(j a (q= O j t); |3 (0,0)>. 



(A2) 

L is the system size, d the dimensionality, j is the para- 
magnetic current operator, and f2„ = 2nnk B T are the 
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bosonic Matsubara frequencies. The diamagnetic contri- 
bution is real and diagonal, and is commonly expressed in 
terms of the carrier density n and the carrier mass m, as 
XapW = ~ S a g(n/m). In interacting Fermi systems, the 
current susceptibilities can be formally represented by two 
classes of diagrams. The first class describes the propaga- 
tion of uncorrelated particle-hole pairs, and can be summed 
to give the so-called "particle-hole bubble", formulated in 
terms of the single-particle spectral function A(k,e). The 
expression of the bubble is 



L ka 



de 1 de 2 Mk,e 1 )A{k,e 2 ) 

x . 

hco + i0 + + e 1 — e 2 



(A3) 



f(e) is the Fermi function, and v k = (l/ft)V£ fc is the group 
velocity with E k the non-interacting electron dispersion. 
The second class of diagrams contains all processes involv- 
ing interactions between the particle and the hole, the so- 
called "vertex corrections", which have been shown to van- 
ish by symmetry in local Fermi liquids characterized by a 
momentum- independent self-energy. 11 This represents a 
considerable simplification, since the exact current-current 
correlation function reduces to Eq. (A3). For a local self- 
energy T,(e) = + iS 2 (e), the spectral function is 



A(k, e y. 



-E 2 0)/7T 



[e-Sfc-SrOO^ + tl^e)] 2 ' 



(A4) 



with £ fc = E k — ix, jLi the chemical potential. Inserting 
Eq. (A4) into Eq. (A3), taking the real part in Eq. (Al), and 
defining the isotropic conductivity as an average over the 

spatial coordinates, a(co) = (l/d)2a=i cr aa( t °)> directly 
leads to Eqs. (8). 

We note in passing that Eqs. (8) provide an interesting 
generalization of the Drude formula for the dc conductivity: 
crjCO) can be transformed into 



de [-fie)] 



2|S 2 (e)| 



du$(e - SxCe) - n|E 2 (e)|) sj^y , (A5) 



(u 2 + l) 2 ' 



with /' the derivative of the Fermi function. Since <3? gen- 
erally is a slow function of its argument, we may expand 
about u = 0, and get at lowest order: 



o-i(O) . 



^ [ - /,W] 2^H e - Si( 4 ^ 



If <& is taken constant, we recover the well-know fact that the 
dc conductivity does not depend on the real part of the self- 
energy explaining the absence of dynamical effective mass 
renormalization in the Drude formula. A weak dependence 
on the dynamical effective mass only enters through the 
energy dependence of the function 



The function in Eq. (8b) is generally a slowly-varying 
function of its argument. After replacing by $(0) in 
Eq. (8a), the £ integration can be performed by means of 
the identity 



dl 



Tir 2 



[( ei -?) 2 + r 2 ][(£ 2 -o 2 + r 2 ] 
ri + r 2 



(e 1 -e 2 ) 2 + [r 1 + r 2 f 



(A7) 



leading to 



*(0) 



de[f{e)-f{e + ha>)} 
-£ 2 0?)-E 2 0? + ftw) 



[hco + Ej(e) - E x (e + hco)] 2 + [E 2 (e) + E 2 (e + hco)] 2 



(A8) 



Equivalently this can be rewritten as: 



a^co) 



HO) r 



de[f(e)-f(e + hco)] 
x Re 



(A9) 



Finally, the requirements that the imaginary part of the con- 
ductivity is positive for positive co, and analytic for Im co > 0, 
lead to Eq. (9) . It can also be seen that the correction to 
Eq. (9), due to the first-order variation of <£(£;), vanishes 
identically if the self-energy has the symmetry property 
E(— e) = — E*(e), as it is the case for the model Eq. (11). 



Appendix B: Calculation of the 5" function 

In order to evaluate the y function, we use the exact 
representation 

1 1 



+00 



-i -rou 

- E 



ip n + x - u ip n — x -u 



(Bl) 



with p n = 2n + 1. Changing variable in the integral in 
(13), we see that each term with p n < is equal to the 
corresponding term with p_ n _ 1 = —p n > 0, so that the sum 
can be rewritten as a sum on n ^ 0. The integrations are 
elementary, 



du 



ip n + x — u ip n — x — u J 1 + x 2 — iy + u 2 
i I 1 1 



r(x,y) yp n + r(x,y) + ix p n + r(x,y)-ix 
r{x,y) = y/l + x 2 -iy, (B2) 
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FIG. 13. Real part (top) and imaginary part (bottom) of the optical 
conductivity normalized by the dc value, as a function of frequency 
and temperature. T is defined in Eq. (11), and a> = 2nk B T a /fi. 
The straight orange lines indicate Hco = 2nk B T, and the green 
lines corresponds to the dome in Fig. 3. 

and we may recast the sums in terms of the digamma func- 
tion ip(z) = liniM^ [logM - J^ =0 l/(n + z)], to finally 
get Eq. (13c). 

Equation (14) shows that the conductivity normalized 

to the dc value is given by (77^, {§jff ) " This is dis " 
played as a function of frequency and temperature in Fig. 13. 
The function 5^ is well approximated by the simple form 
S^i. x ,y) K 1/(1 + x 2 — 3iy/4) over the whole (x,y) plane. 



Appendix C: Comparing Fermi-liquid and Drude formula 

In the Drude regime, defined by hco <C 27ik B T and 
wTqp ~ 1, the conductivity of a Fermi liquid is a uni- 
versal function of given by Eq. (15). This func- 
tion is displayed in Fig. 14. The imaginary part cr 2 i n ~ 
creases first linearly with a slope 1/2 + 3£(3)/n: 2 0.865, 
reaches a maximum at cjt^ fa 1.162, and then decreases as 
12/7i 2 (a)T qp ) _1 . The real part decreases as 16/7i 2 (<uT qp ) -2 , 
and crosses the imaginary part at a>T qp 1.198, slightly 
above the point where o 2 has its maximum. In the Drude 
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FIG. 14. Comparison of the Fermi-liquid conductivity in the 
Drude regime, Eq. (15), with the Drude model, Eq. (1). The 
solid blue and red lines show the real and imaginary parts of 
Eq. (15), respectively. The dashed lines in (a) show the Drude 
model assuming t d = T qp . The dotted lines in (b) show the Drude 
model assuming t d = (n: 2 /12)T qp (see text). 

model Eq. (1), a 2 increases first like cot d , reaches a max- 
imum at cot d = 1, and then decreases as (cdt^) -1 ; <j 1 
decreases as (wt d )" 2 , and crosses ct 2 at its maximum. 

When approximating the low-frequency response of a 
Fermi liquid by the Drude model, it is not possible to choose 
the time t d such that the four remarkable features — initial 
slope and maximum of a 2 , asymptotic decay of <j 1 and 
cr 2 — are all described exactly. If one takes t d = T f , none 
of these features is correctly reproduced, as illustrated in 
Fig. 14a. Matching the initial slope of <j 2 i n b otn mod- 
els implies t d = 0.865T qp , matching the maximum of a 2 
means t d = T qp /1.162, matching the asymptotic decays 
of u l and a 2 leads respectively to t d = (7i/4)Tqp and 
t d = ( 7t2 /12)t (JP . All four determinations of T D /Tqp are 
close to 0.8. The choice t d = (7i 2 /12)T qp , because it gives 
the correct asymptotic decay of u 2 , also ensures that the 
spectral weight is the same in both models by the Kramers- 
Kronig relations. It is furthermore the closest to the average 
of the four determinations. Figure 14b shows a comparison 
of the Fermi-liquid model with the Drude model, assuming 

t d = (7r 2 /12)T 



Appendix D: Memory function in the thermal and Drude 
regimes 

By solving Eq. (2) for M(w), identifying cr(w) with 
the Fermi-liquid result Eqs. (13), and expanding for large 
y = wTqp with the help of Eq. (17), one arrives at the 
expression (21) for the memory function in the thermal 
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regime. Inserting (21) back into Eq. (2), one recovers 
the generalized Drude form Eqs. (18) and (19). From an 
experimental perspective, one sees in (21) that a plot of 
wM 2 ((i))/[w + MjCo))] in the region hco ~ 27ifc B T should 
allow to determine the parameter T . 

We note that, while the real part of M depends on 
the total spectral weight e Q co 2 , the normalized imagi- 
nary part M 2 (w)/M 2 (0) does not, since M 2 (w)/M 2 (0) = 
Recr dc /<7(a>) as shown e.g. in Eq. (7). Unlike M 1 (co), 
M 2 (o>)/M 2 (0) can therefore be determined accurately using 
a low-energy model. With the Fermi-liquid model (13), we 
have 



M 2 (&>) 
M 2 (0) 



Re 



(Dl) 



As lmiy^Retl/^Cx.jO] = (7i/3) 2 (l + x 2 ), one sees that 
the scaling function for M 2 (co) in the thermal regime extrap- 
olates to the value (ti/3) 2 M 2 (0) at co = 0, as can be seen in 
Fig. 5. 

The behavior of M 2 (a>)/M 2 (0) in the Drude regime 
coTqp < 1 is quadratic and given by 



M 2 {co) 
M 2 (0) 



i+ (coz^y 



a + b 



(D2) 



where a and b are the numerical constants a = 1/8 + 
7i 2 /80 - 3a3)/(47r 2 ) - (3C(3)/tt 2 ) 2 m 0.0235 and b = 
1/2 + 7i 2 /60 + 3£(3)/tt 2 w 1.03. This is plotted in Fig. 5a 
as the thin line. 



in the denominator of the integrand, 1 is replaced by 1 +- 
2ZTqpr /ft. The correction term ~ F/T 2 is the ratio between 
the impurity and the electron-electron scattering rates. This 
change can be absorbed into a redefinition of the variable 

y- 



co co + izr/(7ifc B r ) 
y = ^^yi + iy 2 = T2 • 



(El) 



The same modification made in the analytical form Eq. (13c) 
yields the scaling function in the presence of impurity scat- 
tering. It should be noted that, with this modification, the 
function 5" is not unity in the dc limit, and therefore the 
dc conductivity is different from a dc as defined in Eq. (13). 
In particular, the residual dc resistivity can be evaluated as 
p = 2r/[ft*(0)]. 

If co » 2Zr/ft, we can expand in the thermal regime like 
in Eq. (17), which must then be replaced by 



12f 16 

^(x,yi -> oo, y 2 ) = — — +- 



, i l + ^ + -y 2 . (E2) 
71 Ji n y{ V 4 



Consequently, we find that the optical scattering rate reflects 
the increased quasiparticle scattering rate (with an extra 
factor of two), but that the frequency/temperature scaling 
is unchanged: 

/] [[hoof + (27rfc B r) 2 l +■ 2Zr. (E3) 



T opt (w) 3? 



Appendix E: Impurity scattering, robustness of the 2n factor 

In this appendix, we show that the inclusion of a constant 
scattering rate — ir in the Fermi-liquid self-energy Eq. (11) 
does not change the scaling of the optical conductivity in 
the thermal regime. A frequency-independent scattering 
rate can crudely describe the effect of impurity scattering. 
On adding — ir to the self-energy the only change induced 
in Eqs. (13) occurs in the expression of the function 5?: 



The impurity scattering, however, reduces the fre- 
quency/temperature domain below the dome in Fig. 3, 
where T opt (w) > 1/co. The equation of the dome in the 
presence of impurity scattering is 



3 r zr 

Ti (to) = Tr, \ / — I co 



(E4) 



Hence the dome disappears — and consequently the shoulder 
in the log-log plot of cr^co) — if Zr > 3n;fc B r /16. 
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